Dark Matter through the Axion Portal 
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Motivated by the galactic positron excess seen by PAMELA and ATIC/PPB-BETS, we propose 
that dark matter is a TeV-scale particle that annihilates into a pseudoscalar "axion." The positron 
excess and the absence of an anti-proton or gamma ray excess constrain the axion mass and branching 
ratios. In the simplest realization, the axion is associated with a Peccei-Quinn symmetry, in which 
case it has a mass around 360 - 800 MeV and decays into muons. We present a simple and predictive 
supersymmetric model implementing this scenario, where both the Ifiggsino and dark matter obtain 
masses from the same source of TeV-scale spontaneous symmetry breaking. 



Introduction. Evidence for dark matter (DM) is 
by now overwhelming While the precise nature 

and origin of DM is unknown, thermal freezeout of a 
weakly interacting massive particle (WIMP) is a success- 
ful paradigm that arises in many theories beyond the 
standard model. If this is correct, then specific DM prop- 
erties can be probed through direct and indirect detection 
experiments, and pieces of the dark sector might even be 
produced at the Large Hadron Collider (LHC). 

Recent indirect detection results may offer important 
insights into the dark sector. The latest PAMELA 
data 0] is strongly suggestive of a new source of galac- 
tic positrons, bolstering the HEAT [i| and AMS-01 [4] 
anomalies. An intriguing interpretation of the PAMELA 
excess is DM annihilation [1, 0, 0, S]. The rate and 
energy spectrum of the PAMELA positrons are consis- 
tent [7| with the electron/positron excesses seen in the 
balloon experiments ATIC ^ and PPB-BETS [l^, and 
the spectral cutoff in these experiments points to a DM 
mass in the TeV range. 

There are, however, two puzzling features in the 
PAMELA data. First, the positron excess is not ac- 
companied by an anti-proton excess, which strongly con- 
strains the hadronic annihilation modes of the DM 7]. 
Unless DM is as heavy as 10 TeV, the PAMELA data 
disfavors DM annihilation into quarks, Ws, Zs, or Hig- 
gses. Second, the required annihilation cross section in 
the galactic halo is orders of magnitude larger than the 
thermal relic expectation. Therefore, any DM interpreta- 
tion of the data must explain both the large annihilation 
rate and the large annihilation fraction into leptons. 

In this letter, we propose that DM is a TeV-scale parti- 
cle that dominantly annihilates into a pseudoscalar "ax- 
ion" a. The axion mass lies above the electron or muon 
threshold, so that a dominantly decays into leptons with 
suppressed photonic and hadronic decay modes. In this 
scenario, the electron/muon decay channels would ac- 
count for the PAMELA excess, and the photon/pion/tau 
decay channels would be constrained by gamma ray tele- 
scopes like HESS [lH and FERMI [13]. 

In a typical realization of the scenario, DM is a fermion. 
In this case, the dominant annihilation channel is not to 



2a, but to a and a real scalar s. If the scalar domi- 
nantly decays as s —^ aa, each DM annihilation will con- 
tain three axions, yielding a distinctive semi-hard galac- 
tic positron spectrum. The existence of a light scalar s is 
also crucial to enhance the DM galactic annihilation rate 
through nonperturbative effects. 

This DM scenario can arise in any theory where the 
DM mass is generated from spontaneous breaking of a 
global U{l)x symmetry under which leptons have ax- 
ial charges. In two Higgs doublet models, it is natural 
to identify t/(l)x with a Peccei-Quinn (PQ) symmetry 
U{1)PQ [13|, in which case our axion is a heavier variant 



of the DFSZ axion [14[. In this realization, the axion 



must decay dominantly into muons to evade constraints 
from low energy and astrophysical experiments. 

The identification of U{l)x with t/(l)pQ is particu- 
larly well-motived in supersymmetric (SUSY) theories, 
since then both the Higgsino and DM can obtain masses 
from the same source of TeV-scale spontaneous C/(1)pq 
breaking. As we will see, the resulting SUSY model is 
extremely simple and predictive, and leads to interesting 
phenomenology at the LHC. We thus mainly focus on this 
model, although other possibilities are also discussed. 

Basic Setup. To understand the main features of 
the scenario, we first isolate the fields responsible for the 
dominant DM phenomenology and study their dynamics. 
A complete SUSY model will be described later. 

Our starting point is a global U{l)x symmetry that 
is broken by the vacuum expectation value (vev) of a 
complex scalar S: 



(1) 



where a is the "axion," s is a real scalar, and fa is the 
axion decay constant. A vector- like fermion DM ip/^'^ 
obtains a mass from 



C = -^S'V'V' + h.c. 



m-DM — £,fa, 



(2) 



where tp/ip'^ is a standard model (SM) singlet. The sta- 
bility of DM can be ensured by a vector-like symm etry 
acting on i'/ip'^, which could be a remnant of U{l)x [3l| . 
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In order for a to decay into leptons, SM leptons must 
have nontrivial U{l)x charges. For example, in a two 
Higgs doublet model, a coupling of the form 



f{S)hJid + h.c., 



(3) 



can force h^hd, and hence the quarks and leptons, to 
carry nontrivial U{l)x charges. The U{l)x is then a PQ 
symmetry, and we focus on this case until the end of this 
letter. If the coupling of Eq. ^ is sufficiently small, it 
does not drastically affect the DM phenomenology. 

Unlike the ordinary axion [l5| . the mass of a cannot 
come only from pion mixing, since it would then be too 
light to decay into leptons. We therefore need some ex- 
plicit breaking of U{l)x- Also, the potential that gen- 
erates the S vev will also give a mass to s. Both effects 
can be described phenomenologically by the mass terms 
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(4) 



and we assume the hierarchy nia ^ TOs ^ wdm- 
This condition is naturally satisfied in the explicit SUSY 
model considered later. 

The scalar field s decays into aa through the operator 



1 



V2fa 



,s(9^a)2. 



(5) 



arising from the S kinetic term. This is typically the 
dominant decay channel for s. 

Thermal Preezeout. In the above setup, the DM 
has three major annihilation modes 



ss. 



aa, 



sa. 



(6) 



The first two modes do not have an s-wave channel and 
are suppressed in the v Q limit. Annihilation at ther- 
mal freezeout is therefore dominated by the third mode. 

In the limit ms,ma <C m-DM, the w — > annihilation 
cross section is 



[av 



tptp — >sa 



647r/4 



(7) 



A standard thermal relic abundance calculation implies 
1 , , 



[av 



3 X lO^^^'cm-'/s, 



(8) 



so once todm is constrained by future ATIC data, then 
fa is completely determined. As a fiducial value, todm ~ 
1 TeV implies /„ - 1 TeV, and so ^ 1- 

Halo Annihilation. The cross section of Eq. ([S]) is 
too small to account for the observed PAMELA excess. 
For rriDM ^ 1 TeV, the required boost factor is [^, 0] 



(o-w)pAMELA ^ 10 (o-w), 



(9) 



although the precise value is subject to a factor of a few 
uncertainty. Such a large boost factor is difficult to ex- 
plain astrophysically [iq . 



However, the halo annihilation rate can be enhanced 
by nonperturbative effects associated with the light state 
s, with nis <^ TOdm- The relevant effects are the Som- 
mcrfcld enhancement 0, H, [l7| and the formation of DM 
boundstates (WIMPoniums) [18'|. These boost the signal 
by 



B~c , a5 = — , 

nis in 



(10) 



where c is a coefficient which can be as large as 
(ws/TODMi'haio)^ if "^s takes values that allow (near) 
zero-energy boundstates. Here, Whaio ^ 10~^. As we 
will see, our explicit model has nis ~ 0(1 - 10 GeV). 
Combined with a moderate astrophysical boost factor, 
Eq. pH)) can then account for Eq. ^ . The condition for 
the effects being operative in the halo, but not at freeze- 
out, is ?;haio ^ ffroczcout, where Vfreozoout ~ 0.4 [iil. 

The DM annihilation (or para-WIMPonium decay) 
product is sa. After the decay s —^ aa, this yields three 
axions per DM annihilation. There is also an annihila- 
tion channel into tt through s-channel a exchange, but 
its branching fraction is only of 0(1%). This level of 
hadronic activity is consistent with the PAMELA data. 

Since DM does not annihilate directly into leptons, 
the positron injection spectrum is different from [7| and 
closer to [2(2]. If a DM annihilation product comes from 
a (scalar) cascade decay with large mass hierarchies in- 
volving n steps, then its energy spectrum is proportional 
to {ln(mDM/-E')}"~^, where n = for direct annihilation. 
The lepton spectrum per DM annihilation is then 



dE TOnM V 



todm 



E 



1 (£;<mDM). (11) 



This semi-hard lepton spectrum is an unambiguous pre- 
diction of our setup. Since our axion will primarily decay 
into muons, Eq. (jlip must be convoluted with muon de- 
cays to give the final positron spectrum. 

Axion Decays. To account for the observed positron 
excess, the axion must dominantly decay into leptons. 
As we will see, there are strong bounds on the photon 
flux from DM annihilations and this constrains a — > 77 
to be less than k, 1%. Since tt" 77, the decay into 
neutral pions must also be suppressed to the 5% level. 
This disfavors the possibility a t~^t~ , since tan decay 
leads to an 0(1) fraction of 7r"s. 

Compared to the QCD axion, we have an extra 
parameter which affects axion-pion mixing. In terms of 
the mixing angle for the QCD axion O^TrO fn/ fa, the 
new mixing angle is 0aTrO — OaTr'>'m,'^/{m'^ — m?^)- Close 
to the TT^ threshold there is resonant enhancement, but 
for 3> m^, there is a toJ/to^ suppression. A similar 
enhancement also arises for ~ m^. 

For a generic axion, the partial widths to leptons and 
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photons are 



r(a 77) 



16^ ' 
«EM rnl 
1287r3 /2 



(12) 
(13) 



In our case, ce = sin^/3, where tan/3 = {hu)/{hd), and 
depends on vria through 9a-n<^ but is typically 0(1)- 
Except when nia is very far from a lepton mass threshold, 
the photon branching fraction is less than 0(1%). 

The partial width to pions is more complicated. Direct 
decays a — > tttt are suppressed by CP invariance, and 
radiative decays a tttt^ are suppressed by Q!EM/47r. 
The first dangerous channel is a tttttt which arises 
from axion-pion mixing. We estimate the partial width 
as 



r(a 



1 



1287r3 /4 



fir ml 
fa ml 



(14) 



where the combination in parentheses is the approximate 



axion-pion mixing angle when 3> 



Parametri- 



cally, the a — > tttttt mode is two orders of magnitude 
suppressed compared to the a —> mode. As nia 

approaches the yOTT threshold, a — > tttttt is enhanced by 
m?plV^p. Also, a rjiriT decay becomes important around 
the same mass scale, so we estimate the total axion to tt" 
branching ratio to be safe for ma ^ 800 MeV. 

Axion Constraints. The bounds on heavy axions 
are different from ordinary axions. For the range 2mf, < 



< 2m p, a beam-dump experiment at CERN [21 1 



looked for the decay a — > e~'"e~, and definitively rules 
out axion decay constants up to /j, ~ 10 TeV. 

In the region 2m ^ < ma < mx — mj^, our axion decays 
into fJ.'^ n~ with cTa ~ 0(1 ~ 10 /im), and measurements 
of rare kaon decays K — > irfj,'^ fi~ constrain the branch- 
ing ratio K T:a. The estimated branching ratio is 
Bt{K+ TT+a) ^ 3 X 10-6(1 TeV//a)^ and the 

measured rate Br {K^ — > tt+/x+/^-) ~ 1 x lO"'' [is!] is 
consistent with SM expectations. Therefore, this region 
seems to be excluded for /a ~ 1 TeV, especially consid- 
ering that the dimuon invariant mass distribution would 
be peaked at for the axion decay. 

For mii — m^ < ma ^ 800 MeV, there are interesting 
implications for rare T decays. The predicted rate is 
Br(T -> 7a) ~ 3 X 10-6sin'*/3(l TeV//a)^ while the 
experimental bound is Br(T 7a) <^ few x 10"^ for 
prompt a — > fJ.^ fJ-^ decays [24,]. The region of interest 
/a ^ 1 TeV will be tested in future i3-factory analyses. 

In summary, the allowed region for our axion is 



mx — m^ < ma ^ 800 MeV, 



(15) 



and the dominant decay channel is a ^ fi^/i^. For ax- 
ions as heavy as m^ — m-n-, astrophysical bounds are ir- 
relevant. 



Gamma Ray Bound. As already mentioned, the 
axion typically has a non-zero branching fraction into 
photons (or 7t°s), and there are important bounds from 
gamma ray experiments. Since the photon spectrum is 
semi-hard, the strongest bounds come from atmospheric 
Cherenkov telescopes. The expected photon spectrum 
also overlaps with the energy range of FERMI. 

A model-independent constraint comes from a HESS 
study of the Sagittarius dwarf galaxy [2^. They put 
an upper bound on the integrated gamma ray flux for 
> 250 GeV of < 3.6 x lO-^^ cm-^s^^. From 
a given photon spectrum, this can be translated into a 
bound on the DM annihilation cross section 



(av) < 



4TT$^mg)M 
JAil 



250 GcV 



~dE 



■dE 



(16) 



where J ~ 2.2 x 10^^ GeV^ cm^^ is the Sagittarius hne-of- 
sight-integrated squared DM density assuming an NEW 
profile, and Afi = 2 x 10^^ is the HESS solid angle inte- 
gration region. 

Since the annihilation cross section is known from 
Eq. ([5]), we can translate Eq. (fTB|) into a bound on the 
branching fraction to photons. For direct a — > 77 de- 
cays, the energy spectrum is proportional to Eq. (jlip . 
dN-y/dE ~ Br(a — > -fj) dNg/dE, and using the fiducial 
todm = 1 TeV, we obtain the bound 



Br(a ^ 77) ^ 1% (moM = 1 TeV). 



(17) 



The bound on axion decays into pions can be derived 
similarly. Assuming a — > tt'^tt"'"tt-, we obtain 

Br(a ^ tttttt) < 5% (moM = 1 TeV). (18) 

The decay a — > t+t^ leads to a bound on (av) an order 
of magnitude stronger than ((Tw)pamela |33| . 

Supersymmetric Model. In a SUSY context, it is 
natural to assume that the vector-like DM mass is related 
to the vector-like Higgsino mass ij,h- In fact, the simple 
superpotential 



W = C-S**" + XSHuHd, 



(19) 



together with the soft SUSY breaking terms 

Aoft = -e^c - XAxSHaHd - mlS^S+ ■■■ , (20) 

have all the required ingredients, except for the origin of 
the axion mass, which we leave unspecified (can simply 
be a small kS^ term in the superpotential). Without 
the ^ terms, this model is the PQ-symmetric limit of 
the NMSSM, and is sometimes referred to as PQ-SUSY 
[i^ll^]. The vevs for S, Hu and Hd can be generated in 
a stable vacuum, giving m^Ju/^J■H = C/A. 

For A ^ 1 and |m|| <^ X'^v^^, where vew — 174 GeV, 
the dominant phenomenology is determined essentially 
by five parameters 



{mDM, A, tan/3, m?., 771^}, 



(21) 



4 



rriDM 


A 


tan P 


2 


fa 


/iff 


Ax 






ms 




Br(s ^ //) 


ms 


"13/2 


TS 


ma 


Ta 


(JSI [cm^] 


1000 


0.25 


2.0 


-6.8^ 


1100 


270 


650 


110^ 


530^ 


34 


4 • 10-^1 


/ = 6 : 3% 


5.5 


10 eV 


2 ■ 10~^ 


0.7 


8 ■ 10"^^ 


3 ■ 10-« 


1200 


0.10 


4.0 


-6.3^ 


1200 


120 


430 


-69^ 


440^ 


5.6 


1 ■ 10"^^ 


/ = T : 5% 


1.2 


5 eV 


0.02 


0.4 


1 • 10-1* 


4 ■ 10-« 



TABLE I: Two sample spectra in the SUSY model. rriDM, A, tan/?, m|, m3/2, and rria are inputs, and the rest are outputs. 
All the masses are in GeV (except where indicated), and the lifetimes are in seconds, ctsi is a spin-independent DM-nucleon 
cross section, ms is calculated assuming decoupling gauginos. 



with rog and affecting only scalar mixing and axion 
decay, respectively. All other parameters are either de- 
termined by thermal relic calculations, electroweak sym- 
metry breaking, or are secondary to the phenomenology 
relevant here. The model is thus extremely predictive in 
this region, and we present sample spectra in Table HI 

The present SUSY model introduces important addi- 
tions to the minimal structure described before. First, 
the mass of s is no longer a free parameter and is fixed 
by ms ~ AwEwsm(2/3). Second, we have an additional 
light state s, the fermion component of S, whose mass is 
m-s ~ ©(A^wIw/'^SUSy), where msusY is t^H or a gaug- 
ino mass. 

The existence of S states lighter than the electroweak 
scale is consistent with the experimental data. These 
states mix with H^.d states with mixing angles of 
0{v-Ew / fa), and constraints from LEP are satisfied for 
fa Z 1 TeV. Considering tin = Ma AAsin(2/3)/2 
from potential minimization), fa ^ ^ TeV implies A ~ 
0(0. 1), which satisfies the bound on charginos. 

Small values of A allow w 0(1 - 10 GeV), as needed 
for the halo annihilation enhancement. Small A also 
ensures that DM annihilates mainly into S states and 
not Hu.d states, which would give more hadronic activ- 
ity than is allowed by PAMELA. To suppress additional 
hadronic/photonic activity from s decays, the branching 
fraction of s into quarks and taus can be made smaller 
than 0(10%). The s —^ ss mode is subdominant, and 
annihilations of DM into ss are velocity suppressed. 

Since s is hght, it mediates a large DM-nucleon cross 
section, leading to tension with the direct detection 
bound [2^. There are two ways this bound can be satis- 
fied. One is to take to be a few tens of GeV. In this 
case the halo annihilation enhancement occurs through 
(near) zero-energy boundstates. The other is to sup- 
press the s-hd mixing by taking appropriate values of 
|m|| « 0.1A^u|^. In this case the s coupling to the 
nucleon can be accidentally small, with a mild tuning 
of 0(10%). The two cases described here correspond to 
the two points in Table [H The Sommerfeld enhancement 
factors for these points arc > 100 8]. 

For the consistency of the DM story, s must not be 
stable. In low-scale SUSY breaking, it is natural to as- 
sume that s decays into the gravitino G. The mass of 
s is typically above ma, in which case the lifetime is 
given by t~_^^^ ~ 96Trmy^Mp^/m^. For a gravitino mass 



™3/2 ^ 0(10 - 100 eV), this is sufficiently short that s 
never dominates the universe. Also, gravitinos this light 
do not cause a cosmological problem [29| . 

The light s, a, and s states have interesting implica- 
tions for LHC phenomenology. For example, the Higgs 
can decay as ft, ^ aa ^ 4/i. Strongly produced SUSY 
particles will typically cascade decay into the light Hig- 
gsino, which will subsequently decay into s, sometimes 
by emitting an a or s. This leads to pairs and quartets 
of coUimated muons with small invariant mass. 

Finally, we note that in general, DM can either be 
the fermion or scalar component of 5'/^''^, depending on 
which is lighter. Scalar DM works similarly to fermion 
DM, except the dominant annihilation modes are now ss 
and aa. Since the scalar/fermion mass splitting controls 
radiative corrections to m|, the 'i'/'^'^ states could be 
nearly degenerate, leading to coannihilation at freezeout. 

Outlook. We have presented a DM scenario that 
naturally explains the PAMELA and ATIC/PPB-BETS 
data. DM is a TeV-scale particle annihilating into an ax- 
ion a, and the halo annihilation rate is enhanced through 
the scalar s. In the simplest realization, a is associ- 
ated with a PQ symmetry, and we have constructed a 
corresponding SUSY model where the Higgsino and DM 
masses have a common origin from [7(1 )pq breaking. 

There are other implementations of our scenario. For 
example, in models of low-scale dynamical SUSY break- 
ing, the sector breaking SUSY typically leads to an 
R axion with the decay constant fa ^ A/47r, where 
A w 0(10 - 100 TeV) is the dynamical scale. This axion 
can serve as our a if the Higgses and DM obtain U(V}r 
breaking masses. The mass of a is ^ A^/Mpi [30|, 
and the scales could work with 0(1) (or loop) factors. 

One could also consider a purely leptonic axion, for 
example, by introducing a separate U{l)x and Higgs 
fields for the lepton sector. In this case a does not 
have hadronic couplings, eliminating the tension with di- 
rect detection experiments and opening the possibility 
for a — > e"^e" decays. 
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DM may be a Majorana fermion tp obtaining a mass from 
C = -^Stp'^/2+h.c., in which case the stability is ensured 
by a Z2 symmetry. 

Alternatively, we can make V/V''^ charged under a new 
I7(1)ne gauge symmetry. As long as the ?7(1)ne gauge 
coupling is small, (;ne thermal freezeout is un- 

affected and the halo boost is i? ~ QNE/i'haio. How- 
ever, ?7(1)ne must be spontaneously broken to avoid 
constraints from big-bang nucleosynthesis, since axion 
decays would otherwise generate a large energy density 
of ?7(1)ne gauge bosons. Nonperturbative enhancements 
through massless gauge bosons are also ruled out by as- 
trophysical observations ,19]. The model presented in the 
text avoids both bounds automatically, since the scalar s 
has a mass ^ 0(1 GeV). 

The bounds from photon flux may be weakened by as- 
suming a different DM density or velocity profile, or due 
to the uncertainty in Eq. ([5}. In this case, the region 
2mT < rria < 2mt may open up, since hadronic ax- 
ion decays are suppressed by 3mc/mJ tan'*/?. The con- 
straint from rare T decays is satisfied, since the experi- 
mental upper bound on Br(T —^ 7a) with a — > r"'"r~ is 
few X lO"'^ Hi]. The region 800 MeV <, nia <, 2 GeV may 
also be allowed. 



